Radiation from relativistic jets in turbulent magnetic fields 
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Abstract. Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we have investigated 
long-term particle acceleration associated with an relativistic electron-positron jet propagating in an unmagnetized ambient 
electron-positron plasma. The simulations have been performed using a much longer simulation system than our previous 
simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. 
Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. The acceleration of ambient 
electrons leads to a maximum ambient electron density three times larger than the original value. Behind the bow shock in the 
jet shock strong electromagnetic fields are generated. These fields may lead to the afterglow emission. We have calculated the 
time evolution of the spectrum from two electrons propagating in a uniform parallel magnetic field to verify the technique. 
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RPIC SIMULATIONS 

Particle-in-cell (PIC) simulations can shed light on the physical mechanism of particle acceleration that occurs in 
the complicated dynamics within relativistic shocks. Recent PIC simulations of relativistic electron-ion and electron- 
positron jets injected into an ambient plasma show that acceleration occurs within the downstream jet |H 0. [US Ht] . 
In general, these simulations have confirmed that relativistic jets excite the Weibel instability. The Weibel instability 
generates current filaments and associated magnetic fields |7[], and accelerates electrons yj, |2|, |3j U, |5[] . 

Pair Jets Injected into Unmagnetized Pair Plasmas using a Large System 

We have performed simulations using a system with (L x ,L y ,L z ) = (4005A, 105A, 105A) and a total of ~ 1 billion 
particles (12 particles / cell / species for the ambient plasma) in the active grid zones [ 8 ] . In the simulations the electron 
skin depth, X ce = c/©pe = 10. OA, where (Op e = (Ane 2 n e /m s ) l l 2 is the electron plasma frequency and the electron Debye 
length A e is half of the grid size. Here the computational domain is six times longer than in our previous simulations 
Jj,|5|]. The electron number density of the jet is 0.676n e , where n e is the ambient electron density and y = 15. The 
electron/positron thermal velocity of the jet is vj th = 0.014c, where c — 1 is the speed of light. Jets are injected in 
a plane across the computational grid at x = 25A in the positive x direction in order to eliminate effects associated 
with the boundary conditions at x = x m ; n . Radiating boundary conditions were used on the planes at x = x m \ B & x max . 
Periodic boundary conditions were used on all transverse boundaries. The ambient and jet electron-positron plasma 
has mass ratio m e - /m e + — 1 . The electron/positron thermal velocity in the ambient plasma is vjj th = 0.05c. 

Figure 1 shows the averaged (in the y — z plane) electron density and electromagnetic field energy along the jet at 
t = 20000)^' and 37500^' . The resulting profiles of jet (red), ambient (blue), and total (black) electron density are 
shown in Fig. la. The ambient electrons are accelerated by the jet electrons and pile up towards the front part of jet. 
At the earlier time the ambient plasma density increases linearly behind the jet front as shown by the dashed blue line 
in Fig. la. At the later time the ambient plasma shows a rapid increase to a plateau behind the jet front, with additional 
increase to a higher plateau farther behind the jet front. The jet density remains approximately constant except near 
the jet front. 




FIGURE 1. The averaged values of electron density (a) and field energy (b) along the x at t = 37500)^ (solid lines) and 2000a> pe 
(dashed lines). Fig. la shows jet electrons (red), ambient electrons (blue), and the total electron density (black). Fig. lb shows 
electric field energy (blue), magnetic field energy (red), and the total field energy (black) divided by the total kinetic energy. 



The Weibel instability remains excited by continuously injected jet particles and the electromagnetic fields are kept 
at a high level, about four times that seen in a previous much shorter grid simulation system (L x = 640A). At the 
earlier simulation time (t = 2000©^') a large electromagnetic oscillating structure is generated and accelerates the 
ambient plasma. As shown in Fig. lb, at the later simulation time the oscillating structure extends up to x/A = 1 100, 
then becomes more uniform and the magnetic field energy becomes larger than the electric field energy. These strong 
electromagnetic fields become very small beyond xj A = 2000 in the shocked ambient region flBl]. 

The acceleration of ambient electrons becomes visible when jet electrons pass about x/A = 500. The maximum 
density of accelerated ambient electrons is attained at t = 17500)^' . The maximum density gradually reaches a plateau 
as seen in Fig. 1 a. The maximum electromagnetic field energy is located at x/A = 700 as shown in Fig. lb. The location 
of this maximum remains in this region at large simulation times. 

New Numerical Method of Calculating Synchrotron and Jitter Emission from Electron 
Trajectories in Self-consistently Generated Magnetic Fields 

Let a particle be at position ro(f) at time / 10,0]. At the same time, we observe the electric field from the particle 
from position r. However, because of the finite propagation velocity of light, we observe the particle at an earlier 
position ro(t ) where it was at the retarded time t =t — 8t = t — R(t )/c. Here R(t ) = |r — ro(t )| is the distance from 
the charge (at the retarded time t ) to the observer. 

After some calculation and simplifying assumptions the total energy W radiated per unit solid angle per unit 
frequency from a charged particle moving with instantaneous velocity j3 under acceleration j3 can be expressed as 
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Here, n = R(t ) / |R(t )| is a unit vector that points from the particle's retarded position towards the observer. The 
choice of unit vector n along the direction of propagation of the jet (hereafter taken to be the x-axis) corresponds to 
head-on emission. For any other choice of n (e.g., 9 = 1/y), off-axis emission is seen by the observer. The observer's 
viewing angle is set by the choice of n (n 2 + n^ + n\= 1). 

In order to calculate radiation from relativistic jets propagating along the x direction Ja] we consider a test case 
which includes a parallel magnetic field (B x ), and jet velocity of Vji,2 = 0.99c. Two electrons are injected with different 
perpendicular velocities (vxi = 0.1c, v±2 — 0.12c). A maximum Lorenz factor of y max = {(1 — (v? 2 + v^ 2 )/c 2 }~ 1//2 = 
13.48 accompanies the larger perpendicular velocity. 

Figure 2 shows electron trajectories in the x — y plane (a: left panel), the radiation (retarded) electric field (red: Vj_i = 
0.12c, blue: v±\ — 0.1c) (b: middle panel), and spectra (right panel) for the case B x = 3.70. The two electrons are prop- 
agating left to right with gyration in the y — z plane (not shown). The gyroradius is about 0.44A for the electron with the 
larger perpendicular velocity. The power spectra were calculated at the point (x,y,z) = (64, 000, 000A, 43. OA, 43. OA). 
The seven curves show the power spectrum at viewing angles of 0° (red), 10° (orange), 20° (yellow), 30° (moss green), 
45° (green), 70° (light blue), and 90° (blue). The higher frequencies become stronger at the 10° viewing angle . The 
critical angle for off-axis radiation Qy = Ymlx f° r tr, i s case i s 13.48°. As shown in this panel, the spectrum at a larger 
viewing angle (> 20°) has smaller amplitude. 
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FIGURE 2. The paths of two electrons moving helically along the x— direction in a homogenous magnetic (B x ) field shown in the 
x — v-plane (a). The two electrons radiate a time dependent electric field. An observer situated at great distance along the n- vector 
sees the retarded electric field from the moving electrons (b). The observed power spectrum at different viewing angles from the 
two electrons (c). Frequency is in units of co~ l . 



Since the jet plasma has a large velocity x-component in the simulation frame, the radiation from the particles 
(electrons and positrons) is heavily beamed along the x-axis (jitter radiation) fi"ill . 

In order to obtain the spectrum of synchrotron (jitter) emission, we consider an ensemble of electrons selected in 
the region where the Weibel instability has grown fully and electrons are accelerated in the generated magnetic fields. 
We will calculate emission from about 20,000 electrons during the sampling time t s = ?2 — 1 1 with Nyquist frequency 
On = 1 /2At where At is the simulation time step and the frequency resolution A© = l/? s . 

Emission obtained with the method described above is self-consistent, and automatically accounts for magnetic 
field structures on the small scales responsible for jitter emission. By performing such calculations for simulations 
with different parameters, we can investigate and compare the quite different regimes of jitter- and synchrotron- 
type emission [11]. The feasibility of this approach has already been demonstrated I9 ul0ll . and its implementation is 
straightforward. Thus, we should be able to address the low frequency GRB spectral index violation of the synchrotron 
spectrum line of death ifllll . 
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